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Abstract. - Thermal behaviour of superconductors with complex order parameter symmetry 
is studied within a weak coupling theory. It is shown numerically, that the thermal nature of the 
different components of complex order parametrs are qualitatively different. Within the complex 
order parameter scenario, the recent experimental observations by Krishna et ai, [Science 277, 
83 (1997)] on magnetothermal conductivity and by J. Ma et al., [Science 267, 862 (1995)] on 
temperature dependent gap anisotropy for high temperature superconductors can have natural 
explanation. 



An important challenge for current research is to reconcile various conflicting results on sym- 
metry of the superconducting order parameter in high T c cuprates by direct determination of 
the coexistence of other components with d x 2_ y 2 symmetry. Apparently, the long standing 
controversy concerning s-wave versus d-wave pairing in cuprates is gradually turning in favor of 
the d x 2_ y 2. For example, the precise measurements of spontaneously generated half-integral 
flux quanta on bicrystal and trycrystal films E| together with the corner SQUID expermi- 
ments [^) suggest that the superconducting gap changes sign on the Fermi surface. Recent 
measurements on magnetic penetration depth j^], the nuclear spin relaxation rate Q, angle 
resolved photoemission data ^| indeed provide evidence for the existence of corresponding 
low energy excitation states. While there exists some experimental evidence for each of these 
effects H 0, [I it is by no means conclusive and an active debate continues. 

Sun et al. |9j found a nonvanishing tunneling current along the c axis between Y BaiCu%0§+ x 
and the conventional superconductor Pb, which cannot exist in a pure tetragonal d-wave 
superconductor. It has been argued by a number of authors that such data can be explained 
by considering an admixture of s wave component due to orthorhombicity in such materials 
[10 . A self-consistent electronic structure calculation for a d x 2_ y 2 and a d x 2 +y 2 + id xy vortex 
[11 reveals that the scanning tunneling spectroscopy data on vortices in YBa 2 Cu 3 OT-s [|| 
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is inconsistent with simple d x 2_ y 2 symmetry, but consistent with d x i_ y i + id xy . Based on 
such recent important experimental findings, it was suggested by Laughlin J12| that there 
should be a tendency for the high T c superconductors to develop a small d xy order parameter 
(OP) on top of the usual d x 2_ y 2. Direct evidences for such suggestion are provided by the 
recent outstanding experimental data on thermal conductivity by Krishna et al., [ p"3| and 
theoretically, by Wheatly et al, and most recently by Ramakrishnan fl^| . We summarize 
below their essential findings fl3|| as it is found in this work that the complex order parameter 
symmetry has natural explanation to this data. 

A series of high resolution measurements on thermal conductivity (re) in the Bi^Sr-iC aCu-2,0% 
by Krishna et al. show that the re at low temperature becomes field-independent above 
a temperature dependent threshold field Hk(T). Below T c (92 K) and in zero field a broad 
anomaly in re was observed that peaks near 65 K and an applied field (H \\ c) suppresses 
the anomaly. This remarkable result indicates a phase transition separating a low-field state 
where the thermal conductivity decreases with increasing field and a high-field one where it 
is insensitive to applied magnetic field. The authors argue that this phase transition is not 
related to the vortex lattice because of the temperature dependence of the field Hk (T) (which 
is roughly proportional to T 2 ) as well as its magnitude. Instead, they suggest a field-induced 
electronic transition leading to a sudden vanishing of the quasi-particle contribution to the heat 
current. Possible scenarios would be the induction of either id xy or is component with d x 2_ y 2 
symmetry with application of a weak field. Such proposition of complex order parameter 
symmetry in cuprate superconductors reminds us of another important data on angle resolve 
phtoemission experiment (ARPES) by J. Ma et al, [H in which a temperature dependent gap 
anisotropy in the oxygen- annealed Bi2Sr2CaCuO$j rX compound was found. The measured 
gaps along both high symmetry directions (r — M i.e, Cu-0 bond direction in real space and 
T—Xi.e, diagonal to Cu-0 bond) are non-zero at lower temperatures and their ratio is strongly 
temperature dependent. The experimental observation [ pj is however, not reproduced by any 
other group. This observed feature cannot be explained within a simple d-wave scenario and 
has been taken as a signature of a two component order parameter, d x 2_ y 2 type close to T c and 
a mixture of both s and d otherwise Jl5| . While the conductivity data in cuprates E3] indicate 
a finite induction of is (or id xy ) component with d x 2_ y 2 due to application of magnetic field, 
the ARPES data [[l4| indicates they are intrinsically so. 

In this note, we work out the phase diagram of a superconductor with complex order 
parameter symmetries such as d x 2_ y 2 + id xy and d x 2_ y 2 + is (for our purpose) comprising the 
amplitudes of different components of order parameters as a function of the relative pairing 
strength in different channels. It is found that the appearence of d xy component hardly affects 
the d x 2_ y 2 component in the case of a d x 2_ y 2 + id xy symmetry, whereas the occurence of 
s-component strongly suppresses the d x 2_ y 2 gap for d x 2_ y 2 + is symmetry. These effects have 
been substantiated by calculating the temperature dependence of different components (e.g, 
A xy , A s and A4 2 _ 2 ). Interestingly enough, the thermal growth of d x 2_ y 2 is locked at the 
onset of the s-wave component in d x 2 _ y 2 + is symmetry whereas no such strong competition 
is found in case of d x 2 _ y 2 +id xy symmetry. These anomalous thermal behaviors are likely to 
have important impacts on different physical properties like the specific heat, Knight shift etc. 
Within this complex order parameter scenario, we calculate electronic thermal conductivity 
and the temperature dependent gap anisotropy (Ap-Ai / Ap-x Vs. T). Resemblence of the 
calculated results with the observed data are remarkable. 

Assuming the applicability of weak coupling theory to high T c cuprates, the superconducting 
gap equation may be written as 
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where the pairing potential Vkk' for a two component order parameter with a separable form 
and the corresponding gap functions are obtained as, 

2 2 

= E v iflfl> & A * = E ( 2 ) 

For a d K 2_ y 2 + symmetry, Vi = V dx2y2 , V 2 = Vd xy , fl = {cosk x - cosfc a ), fl = 
2 sin k x sin fcy and the corresponding component gap functions are Ai = Ad x2 _ 2 and A2 = 
iAj, . Similarly, for a corresponding d x 2_ y 2 +is phase, V% = V s , f% = constant and A2 = iA s . 
Substituting (2) in (1) and separating the real and imaginary parts, gap equations for different 
components may be written as, 

where the quasiparticle energy spectrum of the superconducting state is given by E k = 
v / (e k -fi) 2 + I A fc | 2 , e k = -2t(cosk x ~\~ cos ky) being the normal state band energy. The 
coupled equations (3) are solved numerically selfconsistently together with a number conserving 
equation to fix chemical potential (//), for a given set of parameters. Then the self-consistent 
values of the order parameters are used to calculate thermal conductivity using the formula 
proposed by Bardeen et al. |lq| , long time ago. 

v (^ fc cosg) 2 dfl 

^ TV 1 dE k > [> 

where v k cos 9 is the component of group velocity parallel to — VT, T is the relaxation rate and 
is the Fermi distribution function. It is important to consider the correct energy dependence 
of F to include strong inelastic and impurity scatterings. Phase diagrams of a d x i_ y i + id xy and 
that of a d x 2_ y 2 +is superconductor evaluated at T — 5 K are presented in the figures 1(a) and 
1(b) respectively. The phase diagram for d x 2_ y 2 + id xy superconductors presented in Fig. 1(a) 
comprises the amplitudes of A x 2_ y 2, A xy as a function of the relative pairing strength between 
the two channels, Vd 2 _ 2 /Vd my - When the relative pairing strength is about 0.65 the A da . v 
starts developing and its appearence affects the A^ 2 2 component only a little, only when 
A xy » Ad 2 2 . A similar phase diagram for the d x 2_ y 2 + is superconductors is presented in 
Fig. 1(b). It is clear from Fig. 1(b) that with the onset of the s wave component, the A^ 2 2 
component is strongly suppressed resulting in a pure s-wave phase for V d 2 2 /V s < 0.39. This 
is because of the peculiar momentum dependence of A^ 2 2 (fc) and A^ ais (k) ; the regions of 
the Fermi surface where A^ 2 2 (fe) is maximum Ad xy (k) is minimum and vice versa. This is 
however, not the case for A^ 2 2 (k) and an isotropic A s symmetry. Note, while obtaining 
phase diagrams in figures 1(a) and 1(b), the stability conditions of different phases are also 
checked through free energy comparison. 

The strong competition between the d x 2_ y 2 and s-wave components in the d x 2_ y 2 + is 
phase is more distinctly evident in Fig. 2(b). In Fig. 2(b) temperature dependence of 
different components of the complex order parameter d + is is presented for various values 
of Vd 2 _ 2 /V s (such that both the order parameters coexist). The curves with left arrow sign 
correspond to thermal dependence of the the s-wave component and the rest that for the 
d-w&ve component. It is shown that with the onset of s-wave component the thermal growth 
of the d-wave component is arrested. This scenario is however absent in case of a d + id 
superconductor (cf. Fig. 2(a)). In Fig. 2(a), with the decrease of the ratio Vd 2 _ 2 jVd the 
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d xy component enhances substantially but the magnitude of d x 2_ y 2 component is reduced only 
marginally. To note, from figures 2(a) and 2(b) that the transition temperature (T c = 84 K) 
is always determined by the d x 2_ y 2 component irrespective of d + id or d + is symmetry of 
the OP. Or in other words, we have d x 2_ y 2 symmetry towards higher temperatures close to T c 
and an admixture of either d + id or d + is otherwise. Furthermore, the only parameter that 
has been tuned through out is the ratio Vd 2 2 /Vd xy and the other cut-off parameter Cl c = 0.6 
(measured with respect to the hopping integral t) is kept fixed. 

According to the proposal jl3| (see also 12), application of the magnetic field induces either 
a is or a id xy component to the d x 2_ y 2 symmetry (i.e stronger the field more the induction 
of complex component) and thereby the magnetothermal conductivity is suppressed at lower 
temperatures with magnetic field. In the complex order parameter scenario, we have a device 
to enhance the complex component (to the d x 2_ y 2 symmetry) by reducing the parameter 
Vd 2 _ 2 /^4» or Vd 2 _ a/Vs- Hence, an essentially similar physical effect can be brought in. 
As a word of caution, however, it is not known whether the magnetic field simply reduces the 
relative pairing strength in a complex order parameter symmetry. In figures 3(a) and 3(b) 
thermal conductivities of a pure superconductor with complex order parameter symmtries (d 
+ id or d + is) are presented for various values of the relative pairing strength between different 
channels. By pure superconductor we mean that the impurity bound states are neglected and 
only essentially anisotropic superconducting states which meet the case of resonance scattering 
is considered. 

Detailed theory of thermal conductivity for unconventional (pure) superconductors are 
discussed by many authors |i~7| . It is worth mentioning that both the d + id and d + is 
symmetries arc parity and time reversal symmetry violating states. Also both the symmetries 
correspond to a fully gapped situation similar to an isotropic s wave superconductor. A 
large suppression in the normalized thermal conductivity i ^'^r) ) at lower temperatures with 
lowering of Vd , , /Vd or Vd 2 , IV s is seen in Fig.s 3(a, b). Apart from a low temperature 
anamoly, a broad maxima is found around T = 61 K below T c = 84 K in both the cases 
(cf Fig. 3(a) and Fig. 3(b)). Also there is only a very little change in the conductivity at 
lower temperatures, when the ratio Vd 2 _ 2 lVd wy or Vd 2 _ 2 /V s is lower enough (for example 
cf. Fig 3(a) for V d ^_ y2 /V dxy = 0.446 and 0.435). For higher values of Vd x2 _ y2 /V dxy =0.555 or 
Vd sc 2_ 2 /Vs =0.455 the low temperature anomaly is worth noticing. The id xy or is component 
being small (for those values of Vd 2 _ =0.555 or Vd 2 _ 2/Vs =0.455) exists only upto 

very lower temperatures (cf. Figs. 2(a) and 2(b) also) and hence it corresponds to a situation 
of a pure <i-wave with nodes at higher temperatures to that of a fully gapped d + id or d + 
is superconductor at lower temperatures. This results in power law falling of thermal conduc- 
tivity at higher temperatures to a sudden exponential fall of the same at lower temperatures 
and hence explains the lower temperature anomaly. However, such temperature dependence 
changes from power law to exponential even at higher temperatures as the complex component 
is or id becomes substantial. These results are in complete agreement with the experimental 
findings p"3[ . There is another extra cusp like feature seen in case of d + is superconductors 
(e.g, Vd 2 _ 2 /V s =0.455) at around T = 14 K. This is reflection to the arresting of the growth 
of d wave component with the appearence of s component seen in Fig. 2(b). However, our 
figures 3(a,b) does not show suppression of the broad peak at around 61 K with field which 
is seen in experiment. This is because the applied field in the experiment is temperature 
dependent (roughly proportional to T 2 ). Therefore, it turns out that the suppression in the 
broad peak may be associated with the temperature field coupled with the magnetic field in 
[ p"3| . So in order to see the net effect of a pure magnetic field in thermal conductivity, it would 
be suggestive to repeat the same experiment |13] with magnetic field which is independent of 
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temperature. Finally, the temperature dependent gap anisotropy i.e, the ratio Ar— Af/Ar-Jc &s 
a function of temperature is plotted in the inset Fig. 1(a) for a superconducting gap function 
with d x 2_ y 2 + id xy symmetry. It is seen that the temperature dependent gap anisotropy 
could be as high as 7 for Vd 2 _ X /Vd xy — 0.446 (at a temperature where the d xy component 
becomes vanishingly small). Note, in the d x 2_ y 2 + \d xy scenario, the Ar-M °c A^ 2 _ 2 and 
A r -x cx A dwy (r - M = (0,tt) and T - X = (tt/2,tt/2)) and therefore, the inset Fig. 1(a) 
can be qualitatively estimated from Fig. 2(a). In inset Fig. 1(a), the large value of the gap 
ratio can be obtained only when the d xy component is comparable to the d x 2_ y 2 component. 
Similar feature is also true in case of d x 2_ y 2 +is superconductors (not shown in figure) [p"5f . 
While the inset Fig. 1(a) would naturally account the experimental obseravtion by Ma et al. 
[|4j, according to Laughlin's conjecture (l^] the d xy component cannot be too large. 

In summary, motivated by recent experimental datas concerning superconducting pairing 
symmetry in high temperature superconductors we illustrated in this letter some interesting 
features of the complex order parameter symmetry. Taking examples of d+id and d+is, it is 
shown that the different components of a complex OP interfer with each other very differently. 
Thermal conductivity is calculated in the complex order parameter phase and found to be in 
agreement with experimental observations based on coupled effect of temperature and magnetic 
field. The present work neither resolves pairing mechanism nor establishes pairing symmetry 
in cuprates but, with reference to recent experimental results [||[ [l4j, the results presented in 
this work may have some important bearings to high-T c cuprate superconductors. 
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Figure Captions 

Fig 1. Phase diagram of superconductors with complex order parameter symmetries (a) 
d x 2_ y 2+id xy and (b) d x 2_ y 2+is. The notations Vd x2 _ y2 , Vd xy , V s refer to pairing strengths in 
the respective channels. The inset Fig. 1(a) represents temperature dependent gap anisotropy 
in the d x 2_ y 2+id xy scenario. 

Fig 2. Thermal variations of different components of the superconductong gap amplitudes 
in the complex order parameter (a) d x 2_ y 2+id xy and (b) d x 2_ y 2+is symmetry for various 
values of Vd 2 2 /Vd xv an d Vd 2 2 /^s respectively. Strong interference between different gap 
components in d x 2_ y 2+is symmetry is worth noticing in contrast to that in d x 2_ y 2+id xy phase. 

Fig. 3 Normalised thermal conductivity (in arbritary units) as a function of temperature is 
shown in case of (a) d x 2_ y 2+id xy and (b) d x 2_ y 2+is pairing symmetry. (The notations k Sj „ 
refer to thermal conductivity in the superconducting and normal state respectively). Loss 
of quasi-particle current with the enhancement of the complex component is seen at lower 
temperatures. 
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